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for transmission. 
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QUALITY BASED IMAGE COMPRESSION 
BACKGROUND OF THE INVENTION 

5 I. Field of the Invention 

The present invention relates to image processing. More specifically, the present 
invention relates to a quality based compression scheme for image signals utilizing 
adaptively sized blocks and sub-blocks of encoded discrete cosine transform coefficient 
10 data. 

II. Description of the Related Art 

In the field of transmission and reception of video signals such as are used for 

15 projecting "films" or "movies", various improvements are being made to image 
compression techniques. Many of the current and proposed video systems make use of 
digital encoding techniques. Digital encoding provides a robustness for the 
communications link which resists impairments such as multipath fading and jamming or 
signal interference, each of which could otherwise seriously degrade image quality. 

20 Furthermore, digital techniques facilitate the use signal encryption techniques, which are 
found useful or even necessary for governmental and many newly developing commercial 
broadcast applications. 

High definition video is an area which benefits from improved image compression 
techniques. When first proposed, over-the-air transmission of high definition video (or 

25 even over-wire or fiber-optical transmission) seemed impractical due to excessive 
bandwidth requirements. Typical wireless, or other, transmission systems being 
designed did not readily accommodate enough bandwidth. However, it has been realized 
that compression of digital video signals may be achieved to a level that enables 
transmission using reasonable bandwidths. Such levels of signal compression, coupled 

30 with digital transmission of the signal, may enable a video system to transmit with less 
power and with greater immunity to channel impairments while occupying a more 
desirable and useful bandwidth. 
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2 

Many compression techniques available offer significant levels of compression, 
but result in a degradation of the quality of the video signal. Typically, techniques for 
transferring compressed information require the compressed information to be transferred 
at a constant bit rate. 

5 One compression technique capable of offering significant levels of compression 

while preserving the desired level of quality for video signals utilizes adaptively sized 
blocks and sub-blocks of encoded Discrete Cosine Transform (DCT) coefficient data. 
This technique will hereinafter be referred to as the Adaptive Block Size Discrete Cosine 
Transform (ABSDCT) method. This technique is disclosed in U. S. Patent No. 
10 5,021,891, entitled "Adaptive Block Size Image Compression Method And System/' 
assigned to the assignee of the present invention and incorporated herein by reference. 
DCT techniques are also disclosed in U. S. Patent No. 5,107,345, entitled 'Adaptive 
Block Size Image Compression Method And System, " assigned to the assignee of the 
present invention and incorporated herein by reference. Further, the use of the ABSDCT 
15 technique in combination with a Differential Quadtree Transform technique is discussed 
inU. S. Patent No. 5,452,104, entitled 'Adaptive Block Size Image Compression Method 
And System/* also assigned to the assignee of the present invention "and incorporated 
herein by reference. The systems disclosed in these patents utilize what is referred to as 
' intra-frame" encoding, where each frame of image data is encoded without regard to the 
20 content of any other frame. Using the ABSDCT technique, the achievable data rate may 
be reduced from around 1.5 billion bits per second to approximately 50 million bits per 
second without discernible degradation of the image quality. 

The ABSDCT technique may be used to compress either a black and white or a 
color image or signal representing the image. The color input signal may be in a Y1Q 
25 format, with Y being the luminance, or brightness, sample, and I and Q being the 
chrominance, or color, samples for each 4x4 block of pixels. Other known formats such 
as the YUV, YC b C y or RGB formats may also be used. Because of the low spatial 
sensitivity of the eye to color, most research has shown that a sub-sample of the color 
components by a factor of four in the horizontal and vertical directions is reasonable. 
30 Accordingly, a video signal may be represented by four luminance components and two 
chrominance components. 
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Using ABSDCT, a video signal will generally be segmented into blocks of pixels 
for processing. For each block, the luminance and chrominance components are passed to 
a block interleaver. For example, a 16x16 (pixel) block may be presented to the block 
interleaver, which orders or organizes the image samples within each 16x16 block to 
5 produce blocks and composite sub-blocks of data for discrete cosine transform (DCT) 
analysis. The DCT operator is one method of converting a time and spatial sampled 
signal to a frequency representation of the same signal. By converting to a frequency 
representation, the DCT techniques have been shown to allow for very high levels of 
compression, as quantizers can be designed to take advantage of the frequency 
10 distribution characteristics of an image. In a preferred embodiment, one 16x16 DCT is 
applied to a first ordering, four 8x8 DCTs are applied to a second ordering, 16 4x4 DCTs 
are applied to a third ordering, and 64 2x2 DCTs are applied to a fourth ordering. 

The DCT operation reduces the spatial redundancy inherent in the video source. 
After the DCT is performed, most of the video signal energy tends to be concentrated in 
15 a few DCT coefficients. An additional transform, the Differential Quad-Tree Transform 
(DQT), may be used to reduce the redundancy among the DCT coefficients. 

For the 16x16 block and each sub-block, the DCT coefficient valoes and the DQT 
value (if the DQT is used) are analyzed to determine the number of bits required to 
encode the block or sub-block. Then, the block or the combination of sub-blocks that 
20 requires the least number of bits to encode is chosen to represent the image segment. For 
example, two 8x8 sub-blocks, six 4x4 sub-blocks, and eight 2x2 sub-blocks may be 
chosen to represent the image segment. 

The chosen block or combination of sub-blocks is then properly arranged in order 
into a 16x16 block. The DCT/DQT coefficient values may then undergo frequency 
25 weighting, quantization, and coding (such as variable length coding) in preparation for 
transmission. Although the ABSDCT technique described above performs remarkably 
well, it is computationally intensive. Thus, compact hardware implementation of the 
technique may be difficult. 



BNSDOCID: <WO 01S6298A1J_> 



WO 01/56298 




PCT/US01/02751 



4 

SUMMARY OF THE INVENTION 

Alternative techniques that make hardware implementation more efficient offer 
certain advantages. Some systems utilize adaptively sized blocks and sub-blocks of 

5 Discrete Cosine Transform (DCT) coefficient data. Although portions of DCT based 
systems utilize quality as a compression parameter, other portions of the data are based 
on the encoding rate, as opposed to using a quality-based metric. An example of such a 
coding rate based paramter is the quantization step selection of the contrast based 
adaptive block size image compression algorithm. 

10 The present invention is a quality based system and method of image 

compression that utilizes adaptively sized blocks and sub-blocks of Discrete Cosine 
Transform coefficient data and a quality based quantization scale factor. A block of pixel 
data is input to an encoder. The encoder comprises a block size assignment (BSA) 
element, which segments the input block of pixels for processing. The block size 

15 assignment is based on the variances of the input block and further subdivided blocks. In 
general, areas with larger variances are subdivided into smaller blocks, and areas with 
smaller variances are not be subdivided, provided the block and sub-block mean values 
fall into different predetermined ranges. Thus, first the variance threshold of a block is 
modified from its nominal value depending on its mean value, and then the variance of the 

20 block is compared with a threshold, and if the variance is greater than the threshold, then 
the block is subdivided. 

The block size assignment is provided to a transform element, which transforms 
the pixel data into frequency domain data. The transform is performed only on the block 
and sub-blocks selected through block size assignment. The transform data then 

25 undergoes scaling through quantization and serialization. Quantization of the transform 
data is quantized based on an image quality metric, such as a scale factor that adjusts with 
respect to contrast, coefficient count, rate distortion, density of the block size 
assignments, and/or past scale factors. Zigzag scanning may also be utilized to serialize 
the data to produce a stream of data. The stream of data may be coded by a variable 

30 length coder in preparation of transmission. The encoded data is sent through a 
transmission channel to a decoder, where the pixel data is reconstructed in preparation for 
display. 
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It is a feature and advantage of the invention to provide a quality based image 
compression system. 

It is another feature and advantage of the invention to allow flexible image quality 
control through management of the bit rate on a frame by frame basis. 
5 It is another feature and advantage of the invention to allow flexible image quality 

control through management of the bit rate on a block by block basis. 

It is another feature and advantage of the invention to maintain quality image 
compression and bit rate control of data accompanying activities such as bursts of 
motion. 

10 It is another feature and advantage of the invention to use signal to noise ratio 

parameters to quantify image quality. 

It is another feature and advantage of the invention to utilize a quantization scale 
factor that adjusts with respect to the contrast of the image. 

It is another feature and advantage of the invention to utilize a quantization scale 
15 factor that adjusts with respect to the AC coefficient count of the DCT blocks that 
comprise an image. 

It is another feature and advantage of the invention to utilize a' quantization scale 
factor that adjusts with respect to the distortion and bit rate between frames. 

It is another feature and advantage of the invention to utilize a quantization scale 
20 factor that adjusts with respect to past quantization scale factors. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The features and advantages of the present invention will become more apparent 
from the detailed description set forth below when taken in conjunction with the 
5 drawings in which like reference characters identify correspondingly throughout and 
wherein: 

FIG. 1 is a block diagram of an encoder portion of a quality based image 
processing system incorporating the variance based block size assignment system and 
method of the present invention; 
10 FIG. 2 is a block diagram of a decoder portion of a quality based image processing 

system incorporating the variance based block size assignment system and method of the 
present invention; 

FIG. 3 is a flow diagram illustrating the processing steps involved in variance 
based block size assignment; 
15 FIGS. 4a, 4b, and 4c illustrate an exemplary block size assignment, the 

corresponding quad-tree decomposition, and the corresponding PQR data. 

FIG. 5 is a graph illustrating the relationship of the quantizer scale factor vs. the 
average bits per pixel rate. 

FIG. 6 is a flow chart of a quantizer scale factor based on an AC coefficient count 
20 of DCT blocks. 

FIG. 7a is a flow chart of a quantizer scale factor based on rate distortion. 

FIG. 7b is a continuation of the flow chart of the quantizer scale factor based on 
rate distortion in FIG 7a. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In order to facilitate digital transmission of digital signals and enjoy the 

5 corresponding benefits, it is generally necessary to employ some form of signal 
compression. To achieve high compression in a resulting image, it is also important that 
high quality of the image be maintained. Furthermore, computational efficiency is desired 
for compact hardware implementation, which is important in many applications. 

The present invention provides a quality based system or apparatus and method 

10 of image compression that takes into account both the image quality and computational 
efficiency in performing image compression. Controlling the bit rate of a truly variable 
rate system based on the goal of maintaining a target bit rate inhibits the purpose of 
maintaining a quality image. Instead, the present invention addresses rate control 
strategies based on quality. Quality-based image compression systems may be on a 

15 block level or on a frame level. The block level system generally uses a greater number of 
encoded bits per frame than frame level control because the block level systems use more 
overhead bits for identification of the particular blocks. 

Before one embodiment of the invention is explained in detail, it is to be 
understood that the invention is not limited in its application to the details of the 

20 construction and the arrangement of the components set forth in the following 
description or illustrated in the drawings. The invention is capable of other embodiments 
and are carried out in various ways. Also, it is understood that the phraseology and 
terminology used herein is for purpose of description and should not be regarded as 
limiting. 

25 The image compression of the present invention is based on discrete cosine 

transform (DCT) techniques, such as that disclosed in co-pending U.S. Patent 
Application "Contrast Sensitive Variance Based Adaptive Block Size DCT Image 
Compression", Serial No. 09/436,085 filed on November 8, 1999, assigned to the assignee 
of the present invention and incorporated herein by reference. Generally, an image to be 

30 processed in the digital domain is composed of pixel data divided into an array of non- 
overlapping blocks, NxN in size. A two-dimensional DCT may be performed on each 
block. The two-dimensional DCT is defined by the following relationship: 
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x(m,n) is the pixel at location (m,n) within an NxM block, and 
X(k,l) is the corresponding DCT coefficient. 



Since pixel values are non-negative, the DCT component X(0,0) is always positive 
10 and usually has the most energy. In fact, for typical images, most of the transform 
energy is concentrated around the component X(0,0). This energy compaction property 
is what makes the DCT technique such an attractive compression method. 

The image compression technique utilizes contrast adaptive coding to achieve 
further bit rate reduction. It has been observed that most natural images are made up of 
15 relatively slow varying flat areas, and busy areas such as object boundaries and high- 
contrast texture. Contrast adaptive coding schemes take advantage. of this factor by 
assigning more bits to the busy areas and less bits to the less busy areas. 

As desicribed in the aforementioned pending patent application Serial No. 
09/436,085, contrast adaptive coding is also useful for reducing the blocking effect. The 
20 blocking effect tends to be more perceptible in busy areas of the image. However, it has 
been realized that the blocking effect is reduced when a smaller sized DCT is used. The 
blocking effect becomes virtually invisible when a 2x2 DCT is used, although the bit per 
pixel performance may suffer. 

Further, contrast adaptive methods utilize intraframe coding (spatial processing) 
25 instead of interframe coding (spatio-temporal processing). Interframe coding inherently 
requires multiple frame buffers in addition to more complex processing circuits. In many 
applications, reduced complexity is needed for actual implementation. Intraframe coding 
is .also useful in a situation that can make a spatio-temporal coding scheme break down 
and perform poorly. For example, 24 frame per second movies can fall into this category 
30 since the integration time, due to the mechanical shutter, is relatively short. The short 
integration time allows a higher degree of temporal aliasing. The assumption of frame to 
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frame correlation breaks down for rapid motion as it becomes jerky, lntraframe coding is 
also easier to standardize when both 50 Hz and 60 Hz power line frequencies are 
involved. Television currently transmits signals at either 50 Hz or 60 Hz. The use of an 
intraframe scheme, being a digital approach, can adapt to both 50 Hz and 60 Hz 
5 operation, or even to 24 frame per second movies by trading off frame rate versus spatial 
resolution. 

For image processing purposes, the DCT operation is performed on pixel data 
that is divided into an array of non-overlapping blocks. Note that although block sizes 
are discussed herein as being NxN in size, it is envisioned that various block sizes may be 

1 0 used. For example, a NxM block size may be utilized where both N and M are integers 
with M being either greater than or less than N. Another important aspect is that the 
block is divisible into at least one level of sub-blocks, such as N/zxN//, N/zxN/y, N/zxM//, 
and etc. where i and j are integers. Furthermore, the exemplary block size as discussed 
herein is a 16x16 pixel block with corresponding block and sub-blocks of DCT 

15 coefficients. It is further envisioned that various other integers such as both even or odd 
integer values may be used, e.g. 9x9. 

FIGs. 1 and 2 illustrate an image processing system 100 incorporating the quality 
based compression system of the present invention. The image processing system 100 
comprises an encoder 104 that compresses a received video signal. The compressed 

20 signal is transmitted using a transmission channel 108, and received by a decoder 112. 
The decoder 112 decodes the received signal into image samples, which may then be 
displayed. 

In general, an image is divided into blocks of pixels for processing. A color signal 
may be converted from RGB space to YC,C 2 space using a RGB to YC,C 2 converter 

25 116, where Y is the luminance, or brightness, component, and Ci and C 2 are the 
chrominance, or color, components. Because of the low spatial sensitivity of the eye to 
color, many systems sub-sample the C, and C 2 components by a factor of four in the 
horizontal and vertical directions. However, the sub-sampling is not necessary. A full 
resolution image, known as 4:4.4 format, may be either very useful or necessary in some 

30 applications such as those referred to as covering "digital cinema." Two possible YC|C 2 
representations are, the Y1Q representation and the YUV representation, both of which 
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are well known in the art. It is also possible to employ a variation of the YUV 
representation known as YCbCr. 

In a preferred embodiment, each of the Y, Cb ? and Cr components is processed 
without sub-sampling. Thus, an input of each component of a 16x16 block of pixels is 
5 provided to the encoder 104. For. illustration purposes, the encoder 104 for the Y 
component is illustrated. Similar encoders are used for the Cb, and Cr components. The 
encoder 104 comprises a block size assignment element 120, which performs block size 
assignment in preparation for video compression. The block size assignment element 120 
determines the block decomposition of the 16x16 block based on the perceptual 

10 characteristics of the image in the block. Block size assignment subdivides each 16x16 
block into smaller blocks in a quad-tree fashion depending on the activity within a 16x16 
block. The block size assignment element 120 generates a quad-tree data, called the PQR 
data, whose length can be between 1 and 21 bits. Thus, if block size assignment 
determines that a 16x16 block is to be divided, the R bit of the PQR data is set and is 

15 followed by four additional bits of Q data corresponding to the four divided 8x8 blocks. 
If block size assignment determines that any of the 8x8 blocks is to be subdivided, then 
four additional bits of P data for each 8x8 block subdivided are added. 

Referring now to FIG. 3, a flow diagram showing details of the operation of the 
block size assignment element 120 is provided. The variance of a block is used as a 

20 metric in the decision to subdivide a block. Beginning at step 202, a 16x16 block of pixels 
is read. At step 204, the variance, v76, of the 16x16 block is computed. The variance is 
computed as follows: 



where N=l 6, and x u is the pixel in the z th row, / h column within the NxN block. At step 
206, first the variance threshold T16 is modified to provide a new threshold T'16 if the 
mean value of the block is between two predetermined values, then the block variance is 
compared against the new threshold, T'16. 




25 
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If the variance vJ6 is not greater than the threshold 77 <5, then at step 208, the 
starting address of the 16x16 block is written into temporary storage, and the R bit of the 
PQR data is set to 0 to indicate that the 16x16 block is not subdivided. The algorithm 
then reads the next 16x16 block of pixels. If the variance vl6 is greater than the threshold 
5 77 <5, then at step 210 ; the R bit of the PQR data is set to 1 to indicate that the 16x16 
block is to be subdivided into four 8x8 blocks. 

The four 8x8 blocks, /=1:4, are considered sequentially for further subdivision, as 
shown in step 212. For each 8x8 block, the variance, v8 h is computed, at step 214. At 
step 216, first the variance threshold T8 is modified to provide a new threshold T'8 if the 
10 mean value of the block is between two predetermined values, then the block variance is 
compared to this new threshold. 

If the variance vS, is not greater than the threshold T8, then at step 218, the 
starting address of the 8x8 block is written into temporary storage, and the corresponding 
Q bit, Qi, is set to 0. The next 8x8 block is then processed. If the variance v£, is greater 
15 than the threshold T8, then at step 220, the corresponding Q bit, £?, , is set to 1 to 
indicate that the 8x8 block is to be subdivided into four 4x4 blocks. 

The four 4x4 blocks,/,— 1 :4, are considered sequentially for further subdivision, as 
shown in step 222. For each 4x4 block, the variance, v4 ij9 is computed, at step 224. At 
step 226, first the variance threshold T4 is modified to provide a new threshold T f 4 if the 
20 mean value of the block is between two predetermined values, then the block variance is 
compared to this new threshold. 

If the variance v4 i} is not greater than the threshold T4 9 then at step 228, the 
address of the 4x4 block is written, and the corresponding P bit, P ij9 is set to 0. The next 
4x4 block is then processed. If the variance v4 {j is greater than the threshold T4, then at 
25 step 230, the corresponding P bit, P ij9 is set to 1 to indicate that the 4x4 block is to be 
subdivided into four 2x2 blocks. In addition, the address of the 4 2x2 blocks are written 
into temporary storage. 

The thresholds 77 d, T8, and T4 may be predetermined constants. This is known 
as the hard decision. Alternatively, an adaptive or soft decision may be implemented. 
30 For example, the soft decision varies the thresholds for the variances depending on the 
mean pixel value of the 2Nx2N blocks, where N can be 8, 4, or 2. Thus, functions of the 
mean pixel values, may be used as the thresholds. 
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For purposes of illustration, consider the following example. Let the 
predetermined variance thresholds for the Y component be 50, 1100, and 880 for the 
1 6\ 16, 8x8, and 4x4 blocks, respectively. In other words, 77 6 = 50, T8 = 1 100, and T4 = 
880. Let the range of mean values be 80 and 100. Suppose the computed variance for the 
5 16x16 block is 60. Since 60 is greater than 776, and the mean value 90 is between 80 and 
100, the 16x16 block is subdivided into four 8x8 sub-blocks. Suppose the computed 
variances for the 8x8 blocks are 1180, 935, 980, and 1210. Since two of the 8x8 blocks 
have variances that exceed T8, these two blocks are further subdivided to produce a total 
of eight 4x4 sub-blocks. Finally, suppose the variances of the eight 4x4 blocks are 620, 

10 630, 670, 610, 590, 525, 930, and 690, with corresponding means values 90, 120, 110, 
115. Since the mean value of the first 4x4 block falls in the range (80, 100), its threshold 
will be lowered to T'4=2Q0 which is less than 880. So, this 4x4 block will be subdivided 
as well as the seventh 4x4 block. The resulting block size assignment is shown in FIG. 
4a. The corresponding quad-tree decomposition is shown in FIG. 4b. Additionally, the 

1 5 PQR data generated by this block size assignment is shown in FIG. 4c. 

Note that a similar procedure is used to assign block sizes for the color 
components C| and C 2 . The color components may be decimated horizontally, 
vertically, or both. 

Additionally, note that although block size assignment has been described as a top 
20 down approach, in which the largest block (16x16 in the present example) is evaluated 
first, a bottom up approach may instead be used. The bottom up approach will evaluate 
the smallest blocks (2x2 in the present example) first. 

Referring back to FIG. 1, the PQR data, along with the addresses of the selected 
blocks, are provided to a DCT element 124. The DCT element 124 uses the PQR data to 
25 perform discrete cosine transforms of the appropriate sizes on the selected blocks. Only 
the selected blocks need to undergo DCT processing. 

The image processing system 100 may optionally comprise DQT element 128 for 
reducing the redundancy among the DC coefficients of the DCTs. A DC coefficient is 
encountered at the top left corner of each DCT block. The DC coefficients are, in 
30 general, large compared to the AC coefficients. The discrepancy in sizes makes it 
difficult to design an efficient variable length coder. Accordingly, it is advantageous to 
reduce the redundancy among the DC coefficients. 
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The DQT element 128 performs 2-D DCTs on the DC coefficients, taken 2x2 at a 
time. Starting with 2x2 blocks within 4x4 blocks, a 2-D DCT is performed on the four 
DC coefficients. This 2x2 DCT is called the differential quad-tree transform, or DQT, of 
the four DC coefficients. Next, the DC coefficient of the DQT along with the three 

5 neighboring DC coefficients within an 8x8 block are used to compute the next level DQT. 
Finally, the DC coefficients of the four 8x8 blocks within a 16x16 block are used to 
compute the DQT. Thus, in a 16x16 block, there is one true DC coefficient and the rest 
are AC coefficients corresponding to the DCT and DQT. 

The transform coefficients (both DCT and DQT) are provided to a quantizer for 

10 quantization. In a preferred embodiment, the DCT coefficients are quantized using 
frequency weighting masks (FWMs) and a quantization scale factor. A FWM is a table 
of frequency weights of the same dimensions as the block of input DCT coefficients. 
The frequency weights apply different weights to the different DCT coefficients. The 
weights are designed to emphasize the input samples having frequency content that the 

1 5 human visual or optical system is more sensitive to, and to de-emphasize samples having 
frequency content that the visual or optical system is less sensitive to. The weights may 
also be designed based on factors such as viewing distances, etc. 

The weights are selected based on empirical data. A method for designing the 
weighting masks for 8x8 DCT coefficients is disclosed in ISO/IEC JTC1 CD 10918, 

20 "Digital compression and encoding of continuous-tone still images - part 1 : Requirements 
and guidelines," International Standards Organization, 1994, which is incorporated herein 
by reference. In general, two FWMs are designed, one for the luminance component and 
one for the chrominance components. The FWM tables for block sizes 2x2, 4x4 are 
obtained by decimation and 16x16 by interpolation of that for the 8x8 block. The scale 

25 factor controls the quality and bit rate of the quantized coefficients. 

Thus, each DCT coefficient is quantized according to the relationship: 



DCT q (iJ) = 



8*DCT(i,j) ^ 1 
fwm(i,j)*q 2 



30 where DCT(ij) is the input DCT coefficient, fwm(i j) is the frequency weighting mask, q 
is the scale factor, and DCTq(i j) is the quantized coefficient. Note that depending on the 
sign of the DCT coefficient, the first term inside the braces is rounded up or down. The 
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DQT coefficients are also quantized using a suitable weighting mask. However, multiple 
tables or masks can be used, and applied to each of the Y, Cb, and Cr componants. 

The block of pixel data and frequency weighting masks are then scaled by a 
quantization scale factor element. In a preferred embodiment, there are 32 scale factors 
5 corresponding to average bit rates. Unlike other compression methods such as MPEG2, 
the average bit rate is controlled based on the quality of the processed image, instead of 
target bit rate and buffer status. In a preferred embodiment, the particular quantization 
scale element is preselected by a selector 130. Although, it should be noted that the 
invention does not require that the quantization scale element be preselected. 

10 In one block based embodiment, the quantization scale factor is based on contrast 

(132). It has been found that the human visual system is more sensitive to high contrast 
areas than it is to low contrast areas. The quantization scale factor based on contrast 
uses a predetermined relationship between bit rate and the quantizer scale factor. In an 
embodiment, the BSA function divides a MxN block into sub-blocks based on contrast 

1 5 such that more subdivisions occur in high contrast blocks and fewer subdivisions occur in 
low contrast blocks. The invention quantifies contrast as the ratio of the square root of 
the variance to the average intensity. Specifically, a quantization scale 'factor is chosen 
according to the BSA as given by the following relationship: 



20 



cr," 



1 < / < 4 



where /?, is the number of bits for the type of block, a /~ is the variance of the i th type 
of block, and p t is the percentage of the f h type of block as obtained form the block size 
assignment function. The quantization scale factor is determined using a predetermined 
25 relationship between the bit rate and the quantizer scale factor. 

. The relationship between bit rate and the quantizer scale factor is illustrated in 
FIG. 5. The relationship between the bit rate and the quantizer scale factor is pre- 
determined by evaluating a group of images to determine the average number of bits 
necessary to describe a pixel. As shown, the quantizer scale factor increases as the 
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average number of bits per pixel decreases. Thus, a specific scale factor Q may be chosen 
for a corresponding value of R,. 

In another block based embodiment, a quantizer scale factor is chosen based on 
the coefficient count of the DCT (136). Typically, the coefficient count system and 
5 method requires one bit per block of overhead. A quantizer scale factor method 300 is 
illustrated in FIG. 6. After the block size assignment, a slice of data comprising NxN 
blocks of data is selected (304). In a preferred embodiment, the slice of data comprises 
blocks of 16 rows by 16 columns of pixel data. A DCT is performed on the block of data 
(308), and further sub-blocks may be designated, as described in reference to FIG. 3. A 
10 threshold for each block of data is determined (312). The threshold may be determined at 
any stage of the process. In a preferred embodiment, the threshold is empirically 
predetermined. Then, the number of AC coefficients that exceed in magnitude the 
predetermined threshold is determined by the following relationship (316): 

M =2i ^X(kJ)/\/\X(kJ)\^T} 
*=o/=o.**o/*o 

15 where T is threshold, represents the DCT coefficient at location (*»0 within an 

NxN block, and M \ s the number of AC coefficients within an NxN block of DCT that 
exceed the threshold. A scale factor is then assigned (320) according to the number of AC 
coefficients exceeding the threshold T. 

After the first block of data within a slice is evaluated and assigned a scale factor, 

20 the next block of NxN data within the slice is selected (324). The scale factor is 
incremented (328), and another DCT is performed (332). A mean square error (MSE) is 
determined between the old and new DCTs (336). The scale factor is then decremented 
(340), and another DCT is performed (344). The mean square error (MSE) is again 
determined between the old and new DCTs (348). The two mean square errors are 

25 compared, and the quantizer factor associated with the lower mean square error is 
assigned to the block of data (352). Then, the cycle repeats for each block of data within 
the slice (356). After a quantizer is assigned for each block of data within a slice, the next 
slice of data is selected and the cycle continues (360). 

Another block based embodiment of determining the quantizer scale factor is a 

30 quantizer scaler based on rate distortion (134), and is illustrated in FIG. 7 (400). 
Generally, the rate distortion quantizer operates in a similar manner to the coefficient 
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count quantizer, but requires more overhead bits depending upon how many scale factors 
are scrutinized. After the block size assignment, a slice of data comprising NxN blocks of 
data is selected (404). In a preferred embodiment, the slice of data comprises blocks of 
16 rows by 16 columns of pixel data. A DCT is performed on the block of data (408), 

5 and further sub-blocks may be designated, as described in reference to FIG. 3. A 
threshold for each block of data is again determined (412) at any stage of the process. An 
average distortion D and an average bit per pixel rate R is also determined (416). Then, 
using the relationship as described with respect to the coefficient count method, the 
number of AC coefficients that exceed in magnitude the predetermined threshold (M) is 

10 determined (420). A scale factor is then assigned (424) according to the number of AC 
coefficients exceeding the threshold T. 

After the first block of data within a slice is evaluated and assigned a scale factor, 
the next block of NxN data within the slice is selected (428). The scale factor is 
incremented (432) from the scale factor of the previous block, and another DCT is 

15 performed (436). A MSE is determined between the old and new DCTs (440). An 
average distortion D and an average bit per pixel rate R corresponding to the quantizer 
scale factor associated with the MSE rate in step 440 is then determined (444). 

The scale factor is then decremented (448) from the scale factor of the previous 
block, and another DCT is performed (452). The mean square error (MSE) is again 

20 determined between the old and new DCTs (456). An average distortion D and an average 
bit per pixel rate R corresponding to the quantizer scale factor associated with the MSE 
rate in step 456 is then determined (460). The difference between the two average 
distortion rates D and the difference between the two bit per pixel rates R are determined 
(464). Then, the two scale factors are determined (468) by taking the ratio of the change 

25 in pixel rate and change in average distortion rate by the following: 



scale factors - 



52X 



The scale factor resulting in the maximum absolute ratio is chosen for the current 
block (472). Then, the cycle repeats for each block of data within the slice (476). After a 
quantizer is assigned for each block of data within a slice, the next slice of data is selected 
30 and the cycle continues (480). 
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The quantizer scale factor may also be selected on a frame by frame basis. 
However, frame based rate control generally achieves lower average bit rates than the 
above described block based rate control. Further, frame based rate control requires more 
use of buffers, as bit transfers occur on a frame-to-frame basis, as opposed to a block-to- 
5 block basis. 

In a frame-by-frame quantizer embodiment, a quantization scale factor may be 
based on the density of various block sized assignments (FIG. 1, block 144). The 
quantizer scale factor for an entire frame of an image is determined based on a metric that 
relates to the visual quality of the image. Visual quality is inferred by determining the 

10 number of 16x16 blocks and the number of 8x8 blocks for a given frame. That is, if a 
given frame has more 8x8 blocks than 16x16 blocks, the image is deemed more complex. 
Accordingly, a more complex image generally requires finer quantization in order to 
maintain quality. Conversely, if a given frame has more 16x16 blocks than 8x8 blocks, 
the image is deemed less complex. Accordingly, a lower quantizer scale factor may be 

15 used and still maintain visual quality. Thus, for each frame a BSA algorithm is 
performed, and the scale factor is based on the following: 

L = " 16 .n g *0 

where n 16 represents the number of 16x16 blocks in a frame, n 8 represents the number of 
8x8 blocks in a frame, and L represents the ratio that exceeds a threshold. In a preferred 
20 embodiment, the threshold is empirically predetermined. The scale factor is then chosen 
based on the value of L and the number of 8x8 blocks for a given frame. 

Alternatively, a frame based quantizer scale factor may be based on the quantizer 
scale factor of a previous frame (148). In this embodiment, a current frame is compressed 
using the scale factor of the previous frame, the scale factor corresponding to the scale 

25 factor index of the previous frame + l, and the scale factor corresponding to the scale 
factor index of the previous frame - 1. During the compression of the current frame, the 
actual scale factor is determined from the BSA information as described with respect to 
the frame-by-frame visual quality method discussed above. The compressed output 
corresponding to the scale factor closest to the actual scale factor of the current frame is 

30 determined, and used along with its corresponding scale factor. That is, 
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[<?(.?,,_,+ 1), otherwise 

where represents the index of scale factor for the previous frame, q (s n ^ +1) 

represents the candidate index scale factor one above the previous frame, ?tv, -1) 
represents the candidate index scale factor one below the previous frame, and <7„ 

5 represents the actual scale factor for the current frame. 

Referring back to FIGS. 1 and 2, the quantized coefficients are provided to a 
zigzag scan seriaiizer 152. The serializer 152 scans the blocks of quantized coefficients in 
a zigzag fashion to produce a serialized stream of quantized coefficients. A number of 
different zigzag scanning patterns, as well as patterns other than zigzag may also be 

10 chosen. A preferred technique employs 8x8 block sizes for the zigzag scanning, although 
other sizes may be employed. 

Note that the zigzag scan serializer 152 may be placed either before or after the 
quantizers 132, 136, 140, 144 or 148. The net results are generally the same. 

In any case, the stream of quantized coefficients is provided to a variable length 

15 coder 156. The variable length coder 156 may make use of run-length encoding of zeros 
followed by Huffman encoding. This technique is discussed in detail in aforementioned 
U.S. Pat. Nos. 5,021,891, 5,107,345, and 5,452,104, and is summarized herein. A run- 
length coder would take the quantized coefficients and separate out the zero from the 
non-zero coefficients. The zero values are referred to as run-length values, and are 

20 Huffman encoded. The non-zero values are separately Huffman encoded. 

A modified Huffman coding of the quantized coefficients is also possible and is 
used in the preferred embodiment. Here, after zigzag scanning, a run-length coder will 
determine the run-length/size pairs within each 8x8 block. These run-length/size pairs are 
then Huffman encoded. 

25 Huffman codes are designed from either the measured or theoretical statistics of 

an.image. It has been observed that most natural images are made up of flat or relatively 
slowly varying areas, and busy areas such as object boundaries and high-contrast texture. 
Huffman coders with frequency-domain transforms such as the DCT exploit these 
features by assigning more bits to the busy areas and fewer bits to the flat areas. In 
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general, Huffman coders make use of look-up tables to code the run-length and the non- 
zero values. Multiple tables are generally used, with 3 tables being preferred in the 
present invention, although 1 or 2 can be employed, as desired. 

The compressed image signal generated by the encoder 104 may be temporarily 

5 stored using a buffer 160, and then transmitted to the decoder 112 using the transmission 
channel 108. The PQR data, which contains the block size assignment information, is 
also provided to the decoder 112. The decoder 112 comprises a buffer 164 and a variable 
length decoder 168, which decodes the run-length values and the non-zero values. 

The output of the variable length decoder 168 is provided to an inverse zigzag 

10 scan serializer 172 that orders the coefficients according to the scan scheme employed. 
The inverse zigzag scan serializer 172 receives the PQR data to assist in proper ordering 
of the coefficients into a composite coefficient block. 

The composite block is provided using a selector 174 to an inverse quantizer 176, 
178, 180, 182, or 184, corresponding to the quantizer used in the encoder 104, for 

15 undoing the processing due to the use of the quantizer scale factor and the frequency 
weighting masks. 

The coefficient block is then provided to an IDQT element 186, followed by an 
IDCT element 186, if the Differential Quad-tree transform had been applied. Otherwise, 
the coefficient block is provided directly to the IDCT element 190. The IDQT element 

20 186 and the IDCT element 190 inverse transform the coefficients to produce a block of 
pixel data. The pixel data may then have to be interpolated, converted to RGB form, and 
then stored for future display. 

Accordingly, a system and method is presented for image compression that 
performs block size assignment based on pixel variance. Variance based block size 

25 assignment offers several advantages. Because the Discrete Cosine Transform is 
performed after block sizes are determined, efficient computation is achieved. The 
computationally intensive transform need only be performed on the selected blocks. In 
addition, the block selection process is efficient, as the variance of pixel values is 
mathematically simple to calculate. Still another advantage of variance based block size 

30 assignment is that it is perceptually based. Pixel variance is a measure of the activity in a 
block, and provides indication of the presence of edges, textures, etc. It tends to capture 
the details of a block much better than measures such as the average of pixel values. 
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Thus, the variance based scheme of the present invention assigns smaller blocks to 
regions with more edges and larger blocks to the flatter regions. As a result, outstanding 
quality may be achieved in the reconstructed images. 

Another important advantage is that the block size assignment is made regardless 
5 of bit rate and relative to several quantization scale factors, allowing for greater accuracy 
and clarity in controlling image quality. By allowing for a quality based quantization 
scale factor, image compression is at a truly variable bit rate. This allows the image 
compression system to preserve details in all areas that are above just noticeable 
visibility threshold. Furthermore, the variance based image compression provides a less 
10 noticeable degradation of image quality, unlike methods such as Moving Picture Experts 
Group (MPEG) compression techniques. This is particularly crucial for applications 
such as in the area of digital cinema. 

With digital video in high demand, piracy is a serious threat. Digital watermarking 
is an important requirement to deter copyright infringement, and loss of revenue. As 
15 watermarking is done in areas of an image that are perceptibly significant, variance based 
block size assignment is a natural candidate for use in watermarking. 

The previous description of the preferred embodiments is provided to enable any 
person skilled in the art to make or use the present invention. The various modifications 
to these embodiments will be readily apparent to those skilled in the art, and the generic 
20 principles defined herein may be applied to other embodiments without the use of the 
inventive faculty. Thus, the present invention is not intended to be limited to the 
embodiments shown herein but is to be accorded the widest scope consistent with the 
principles and novel features disclosed herein. 

Other features and advantages of the invention are set forth in the following 

25 claims. 

What we claim is: 
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CLAIMS 



1. A method of compressing a digital image, the image comprising pixel data, 
2 the method comprising the acts of: 

reading a block of pixel data; 
4 generating a block size assignment to divide blocks of pixel into sub-blocks of 

pixel data; 

6 transforming the sub-blocks of pixel data into corresponding frequency domain 

representations; and 

8 scaling the frequency domain representations into a stream of data, wherein the 

act of scaling is based on a quality metric correlating with the quality of the image. 

2. The method of claim 1 , wherein the act of scaling further comprises the act 
2 of providing a frequency weighted mask to said sub-blocks of pixel data, such that the 

frequency weighting mask provides emphasis to the portions of the image that the human 
4 visual system is more sensitive, and provides less emphasis to the portions of the image 
that the human visual system is less sensitive. 

3. The method of claim 1 , wherein the act of scaling further comprises the act 
2 of quantizing the sub-blocks of pixel data based on the quality of the image. 

4. The method of claim 1, wherein the quality metric is the signal to noise 

2 ratio. 

5. The method of claim 1, wherein the act of transforming performs a 
2 Discrete Cosine Transform. 



6. The method of claim 1, wherein the act of transforming performs a 
2 Discrete Cosine Transform followed by a Differential Quad-tree Transform. 



7. The method of claim 5, wherein the act of scaling is on a block by block 

2 basis. 
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8. The method of claim 7, wherein the act of scaling is based on the contrast 
2 of the sub-block of pixel data. 

9. The method of claim 7, wherein the act of scaling is based on the number 
2 of AC coefficients within a sub-block that exceed in magnitude a predetermined 

threshold. 

10. The method of claim 7, wherein the act of scaling is based on the mean 
2 square error of the sub-block. 

1 1 . The method of claim 5, wherein the act of scaling is on a frame by frame 

2 basis. 

12. The method of claim 1 1, wherein the act of scaling is based on the number 
2 and size of block size assignments per frame. 

13. The method of claim 11, wherein the act of scaling is based on the scaling 
2 factor of the previous frame. 

14. The method of claim 1, wherein the act of generating further comprises the 

2 acts of: 

determining a variance of pixel values for said block of pixel data; 
4 comparing said variance with a threshold; 

making a decision to subdivide the block based on a result of the act of comparing; 
6 if said decision is to subdivide said block, then repeating the acts of determining, 

comparing, and making for each sub-block until a predetermined criteria is satisfied; and 
8 designating as the block size assignment each block or sub-block that is not 

further subdivided. 

15. The method of claim 14, wherein the step of making a decision calls for a 
2 block to be subdivided if the variance is greater than the threshold. 

1 6. The method of claim 1 5, wherein the threshold is predetermined. 
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1 7. The method of claim 1 5, wherein the threshold is a function of the mean of 
2 pixel values of the block or sub-block being evaluated. 

18. The method of claim 15. wherein the threshold changes for each level of 
2 subdivision. 

19. The method of claim 14, wherein the predetermined criteria for no longer 
2 repeating the steps of determining, comparing, and making is based on a pre-selected 

minimum block size of pixel data. 

20. The method of claim 1, wherein the act of scaling the coefficient data into 
2 a stream of data further comprises zigzag scanning the coefficient data. 

21. The method of claim 21, wherein the act of zigzag scanning is performed 
2 using an 8x8 pixel block size. 

22. The method of claim 1, wherein the act of scaling further comprises the act 
2 of coding the stream of data into a serialized stream of data in preparation for 

transmission. 

23. The method of claim 22, wherein the act of coding comprises a Huffman 

2 coder. 

24. The method of claim 23, wherein the Huffman coder employs multiple 
2 look-up tables to code run-length and non-zero values. 

25. The method of claim 24, wherein there are three look-up tables. 

26. The method of claim 1, wherein the act of generating is based on the 
2 perceptual characteristics of the image. 
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27. A quality based image compression system for compressing pixel data, the 
2 system comprising: 

block size assignment means to divide blocks of pixel into sub-blocks of pixel 

4 data; 

transform means to transform the sub-blocks of pixel data into corresponding 
6 frequency domain representations; and 

scaling means to scale the frequency domain representations into a stream of data, 
8 wherein the scaling means is based on a quality metric correlating with the quality of the 
image. 

28. The system of claim 27, wherein the scaling means further comprises a 
2 frequency weighted masking means to said sub-blocks of pixel data, such that the 

frequency weighted mask provides emphasis to the portions of the image that the human 
4 visual system is more sensitive, and provides less emphasis to the portions of the image 
that the human visual system is less sensitive. 

29. The system of claim 27, wherein the scaling means further comprises a 
2 quantizer to quantize the sub-blocks of pixel data based on the quality of the image. 

30. The system of claim 27, wherein the quality metric is the signal to noise 

2 ratio. 

31 . The system of claim 27, wherein the transform means performs a Discrete 
2 Cosine Transform. 

32. The system of claim 27, wherein the transform means performs a Discrete 
2 Cosine Transform followed by a Differential Quad-tree Transform 

33. The system of claim 31, wherein the scaling means scales the data on a 
2 block by block basis. 
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34. The system of claim 33, wherein the scaling means is based on the 
2 contrast of the sub-block of pixel data. 

35. The system of claim 33, wherein the scaling means is based on the number 
2 of AC coefficients within a sub-block that exceed in magnitude a predetermined 

threshold. 

36. The system of claim 33, wherein the scaling means is based on the mean 
2 square error of the sub-block. 

37. The system of claim 31, wherein the scaling means is on a frame by frame 

2 basis. 

38. The system of claim 37, wherein the scaling means is based on the number 
2 and size of block size assignments per frame. 

39. The system of claim 37, wherein the scaling means is based on the scaling 
2 factor of the previous frame. 

40. The system of claim 27, wherein the block size assignment means 
2 determines a variance of pixel values for the block of pixel data; compares the variance 

with a threshold; makes a decision to subdivide the block based on a result of the 
4 comparing; wherein if the decision is to subdivide the block, then repeats the 

determination of variance, the comparison with a threshold, and the decision to subdivide 
6 for each subdivided block until a predetermined criteria is satisfied; and designates as the 

block size assignment each block or sub-block that is not further subdivided. 

41 . The system of claim 40, wherein the decision to subdivide calls for a block 
2 to be subdivided if the variance is greater than the threshold. 

42. The system of claim 41, wherein the threshold is predetermined. 
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43. The system of claim 41, wherein the threshold is a function of the mean of 
pixel values of the block or sub-block being evaluated. 



44. The system of claim 41, wherein the threshold changes for each level of 
2 subdivision. 

45. The system of claim 40, wherein the predetermined criteria for no longer 
2 subdividing is based on a pre-selected minimum block size of pixel data. 

46. The system of claim 27, wherein the scaling further comprises a zigzag 
2 scanner. 

47. The system of claim 46, wherein the zigzag scanner employs an 8x8 pixel 
2 block size for the zigzag scanner. 

48. The system of claim 47, wherein the scaling means further comprises a 
2 coder to code the stream of data into a serialized stream of data in preparation for 

transmission. 

49. The system of claim 48, wherein the coder is a Huffman coder. 

50. The system of claim 49, wherein the Huffman coder employs multiple 
2 look-up tables to code run-length and non-zero values. 

51 . The system of claim 50, wherein there are three look-up tables. 

52. The system of claim 27, wherein the generating means is based on the 
2 perceptual characteristics of the image. 

53. A method of compressing a digital image, the image comprising pixel data, 
2 the method comprising the acts of: 

reading a block of pixel data; 
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4 generating a block size assignment to divide blocks of pixel into sub-blocks of 

pixel data; 

6 transforming the sub-blocks of pixel data into corresponding frequency domain 

representations; and 

8 scaling the frequency domain representations into a stream of data, wherein the 

act of scaling is based on a quality metric correlating with the quality of the image; 
10 wherein the act of scaling further comprises the act of providing a frequency weighted 

mask to said sub-blocks of pixel data such that the frequency weighted mask provides 
12 emphasis to the portions of the image that the human visual system is more sensitive, 

and provides less emphasis to the portions of the image that the human visual system is 
14 less sensitive; and wherein the act of scaling further comprises the act of quantizing the 

sub-blocks of pixel data based on the quality of the image. 

54. The method of claim 53, wherein the quality metric is the signal to noise 

2 ratio. 

55. The method of claim 53, wherein the act of transforming performs a 
2 Discrete Cosine Transform. 

56. The method of claim 53, wherein the act of transforming performs a 
2 Discrete Cosine Transform followed by a Differential Quad-tree Transform 

57. The method of claim 53, wherein the act of scaling is on a block by block 

2 basis. 

58. The method of claim 57, wherein the act of scaling is based on the contrast 
2 of the sub-block of pixel data. 

59. The method of claim 57, wherein the act of scaling is based on the number 
2 of AC coefficients within a sub-block that exceed in magnitude a predetermined 

threshold. 
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60. The method of claim 57, wherein the act of scaling is based on the mean 
2 square error of the sub-block. 

61. The method of claim 55, wherein the act of scaling is on a frame by frame 

2 basis. 

62. The method of claim 61, wherein the act of scaling is based on the number 
2 and size of block size assignments per frame. 

63. The method of claim 61 , wherein the act of scaling is based on the scaling 
2 factor of the previous frame. 

64. A quality based image compression system for compressing pixel data, the 
2 system comprising: 

block size assignment means to divide blocks of pixel into sub-blocks of pixel 

4 data; 

transform means to transform the sub-blocks of pixel data into corresponding 
6 frequency domain representations; and 

scaling means to scale the frequency domain representations into a stream of data, 
8 wherein the scaling means is based on a quality metric correlating with the quality of the 
image; wherein the scaling means further comprises a frequency weighted masking means 
10 to said sub-blocks of pixel data, such that the frequency weighted mask provides 
emphasis to the portions of the image that the human visual system is more sensitive, 
12 and provides less emphasis to the portions of the image that the human visual system is 
less sensitive; and wherein the scaling means further comprises a quantizer to quantize 
14 the sub-blocks of pixel data based on the quality of the image. 

64. The system of claim 63, wherein the quality metric is the signal to noise 

2 ratio. 

65. The system of claim 63, wherein the transform means performs a Discrete 
2 Cosine Transform. 
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66. The system of claim 63, wherein the transform means performs a Discrete 
Cosine Transform followed by a Differential Quad-tree Transform 

67. The system of claim 66, wherein the scaling means scales the data on a 
block by block basis. 

68. The system of claim 66, wherein the scaling means is based on the 
2 contrast of the sub-block of pixel data. 

69. The system of claim 66, wherein the scaling means is based on the number 
2 of AC coefficients within a sub-block that exceed in magnitude a predetermined 

threshold. 

70. The system of claim 66, wherein the scaling means is based on the mean 
2 square error of the sub-block. 

71. The system of claim 66, wherein the scaling means is on a frame by frame 

2 basis. 

72. The system of claim 71, wherein the scaling means is based on the number 
2 and size of block size assignments per frame. 
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